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Abstract 


The  upwards  vertical  flow  of  oil-water  and  air-water  mixtures 

has  been  investigated  in  22.88  foot  lengths  of  1.025  inch  and  1.50 

inch  tubing  respectively.  Pressure  drop  data  were  obtained  from 
the  1.025  inch  test  section  up  to  a  maximum  second  phase  (oil)  rate 
of  11.8  ft.  per  sec.  and  a  maximum  constant  superficial  reference 
phase  velocity  of  5*6  ft.  per  sec.  In  the  1.50  inch  tube  second 
phase  densities  of  0.092,  O.362  and  0.552  lbs.  per  cu.  ft.  were 
investigated  with  second  phase  rates  varying  from  1.34  to  46.8  ft.  per 

sec.  and  reference  phase  rates  varying  from  O.0695  to  7*35  ft.  per  sec. 

The  relationship  between  the  input  volume  ratio  of  the  phases  and 
the  pressure  drop  exibited  a  minimum,  a  maximum  and  a  second  minimum 
for  the  air-water  system.  This  also  occurred  above  a  certain  reference 
phase  velocity  for  the  oil-water  system.  Below  a  superficial  reference 
phase  velocity  of  one  foot  per  second  the  maximum  and  the  second  minimum 
were  not  evident  in  the  oil-water  system. 

A  friction  factor,  f^,  based  on  the  superficial  reference  phase 
velocity,  was  calculated  from  equations  developed  by  Govier,  Radford 
and  Dunn  (8)  and  correlated  with  the  second  phase  density,  tube  diameter 
and  velocities.  An  increase  in  second  phase  density  decreases  the  value 
of  f ',  other  variables  remaining  constant. 

■U 

Pressure  drop  regimes,  defined  by  previous  workers  through  the 
loci  of  the  minima  and  maximum  points  of.  the  pressure  drop  curves, 
were  used  to  describe  the  flow  patterns.  Second  phase  density  was  found 
to  affect  the  location  of  all  loci  and  the  transitions  between  flow 


patterns . 
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The  holdup  ratio  was  correlated  with  superficial  reference  phase 
velocity ,  tube  diameter  and  second  phase  density.  The  ratio  was  found 
to  decrease  with  an  increase  in  the  second  phase  density. 
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I  Introduction. 

"Two  phase  flow"  refers  to  the  simultaneous  flow  of  one  or  more 
components  in  two  different  phases.  Usually  the  term  is  applied  to  the 
flow  of  gas-liquid  systems.  It  can  however,,  apply  equally  as  well  to 
liquid' -liquid,  solid-liquid  or  solid-gas  systems.  Problems  in  two  phase 
flow  appear  in  diversified  fields. 

A  few  examples  are: 

(1)  Petroleum  Production 

(a)  Gas  Lift  Installations  (gas-liquid) 

(b)  Plow  through  Porous  Media  (gas-liquid-liquid) 

(c)  Pipeline  Transportation  of  Gas-Oil  Mixtures  (gas-liquid) 

(2)  Process  Industries 

(a)  Evaporators  (vapor-liquid) 

(b)  Refrigerators  (vapor-liquid) 

(c)  Steam  Generators  (vapor -liquid) 

(d)  Catalytic  Cracking  Units  (solid-gas) 

(e)  Air  Lift  Installations  (gas-liquid) 

(3)  Wuclear  Energy 

(a)  Boiling  Water  Reactors  (vapor-liquid) 

(b)  Homogeneous  Aqueous  Reactors  (solid-liquid) 

As  it  appears  presently  impossible  to  treat  all  two  phase  systems 
similarily,  the  discussion  will  be  limited  to  the  gas-liquid  and  liquid- 
liquid  systems  which  exhibit  similar  behaviour. 

In  studies  of  such  two  phase  systems  the  primary  object  is  the 
prediction  of  the  effect  of  changing  the  variables  of  geometry,  flow 
rates,  densities,  viscosities  and  surface  tension  on  the  pressure  drop, 
flow  pattern  and  relative  phase  velocity  or  holdup.  Physical  models  have 


been  developed  by  various  investigators  in  an  effort  to  devise 
suitable  correlation  methods.  In  almost  every  case  the  analogy  to 
single  phase  flow  is  used  as  the  model.  Common  among  correlating 
methods  for  the  pressure  drop  relationships  is  the  friction  factor 
approach.  This  method  can  be  further  subdivided  into  a  friction  factor 
based  on  the  average  properties  of  the  flowing  mixture  and  one  based 
on  one  of  the  phases  alone.  The  latter  approach  has  been  employed 
at  the  University  of  Alberta  in  the  study  of  the  upwards  vertical  flow  of 
air -water  mixtures. 

Studies  of  two  phase  flow  at  the  University  of  Alberta  began  in 
1948  when  Radford  (15)  investigated  the  effect  of  the  air  and  water 
flow  rates  on  the  flow  patterns.  Later ,  Dunn  (4)  continued  the  work 
of  Radford  and  was  able  to  determine  the  effect  of  the  air  and  water 
flow  rates  on  the  pressure  drop,  flow  patterns  and  holdup. -  Results  of 
the  work  of  Radford  and  Dunn  were  published  by  Govier,  Radford  and  Dunn  (8). 
In  1957  Short  (l6)  expanded  the  experimental  facilities  used  by  Radford 
and  Dunn,  and  was  able  to  determine  the  effect  of  tubing  diameter  on  the 
pressure  drop,  flow  patterns  and  holdup  in  two  phase  flow.  Data  collect¬ 
ed  by  Short,  together  with  applicable  data  from  the  work  of  Radford  and 
Dunn,  is  being  published  by  Govier  and  Short  (9)* 

The  present  work  is  an  investigation  of  the  effect  of  the  second 
(or  non-water)  phase  density  . on  pressure  drop,  flow  patterns  and  holdup. 

The  two  phase  system  studied  was  composed  of  water  as  the  reference  phase 
with  air  or  oil  as  the  second  phase .  Three  different  second  phase  den¬ 
sities  were  studied  by  varying  the  average  pressure  in  the  air -water 
system.  A  higher  second  phase  density  was  obtained  with  a  low  viscosity, 
low  density  oil. 

A  study  is  now  underway  on  the  effect  of  the  second  phase  viscosity 
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on  the  pressure  drop,  flow  pattern  and  holdup.  Results  of  this  study 
will  he  used  to  verify  the  effect  of  viscosity  in  the  correlations  of 
flow  pattern  and  pressure  drop  presented  herein. 
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II  Literature  Review. 

Initial  investigators  in  the  field  of  gas-liquid  flow  were  concerned 
mainly  with  the  prediction  of  the  pressure  drop,  although  the  problem  of 
relative  phase  velocity  or  slip  was  also  recognized.  Later,  investigators 
realized  the  importance  of  the  flow  patterns  and  attempted  to  describe 
them  in  relation  to  the  pressure  drop  and  slip.  Recent  studies  of  two 
phase  flow  generally  account  for  the  effect  of  the  basic  variables  on  all 
three  factors;  pressure  drop,  flow  pattern  and  slip  or  "holdup". 

A.  Pressure  Drop.' 

.  The  first  application  of  two  phase  flow  to  receive  attention  was  the 

air-lift  method  of  pumping  fluids.  Davis  and  Weidner  (l)  of  the  University 

of  Wisconsin  were  among  the  first  investigators  to  publish  air-lift  data. 

They  realized  that  there  was  a  necessity  for  performance  data  (pressure  drop 

in  effect)  for  the  air-lift  pump  on  which  accurate  design  could  be  based. 

Experiments  were  conducted  in  a  l-l/4  inch  pipe  with  eduction  pipe  lengths 

varying  from  19  feet  to  42  feet.  Air  rates,  water  rates  and  efficiencies 

were  recorded  as  well  as  the  "submergence  ratio"  or  ratio  of  the  submerged 

pipe  to  total  pipe  length.  Since  the  submergence  ratio  is  equivalent  to 

the  term  vR  AP  in  equation  (6),  in  effect,  pressure  drop  data  were 
AX 

determined.  Reference  is  made  later  to  selected  values  of  Davis  and 
Weidner’ s  data. 

In  1929  Uren  et  al  (19)  began  experimental  studies  of  the  pressure 
drop  associated  with  the  upwards  vertical  flow  of  air  and  an  l8°  A.P.I. 
petroleum  fraction.  Tests  were  conducted  in  a  simulated  gas-lift  well 
consisting  of  42  feet  of  5“l/2  inch  casing  enclosing  4l  feet  of  2  inch 
test  section  tubing.  Air-oil  ratios  varied  from  92.6  to  1250  cubic  feet 
per  cubic  feet  (air  at  standard  conditions).  The  Fanning  friction  factor 
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expression  was  used  to  calculate  a  friction  factor  -which  was  then  plotted 
as  a  function  of  the  Reynolds  number.  The  density  in  the  friction  factor 
calculation  and  the  density  and  viscosity  in  the  Reynolds  number  were 
calculated  using  average  properties  of  the  insitu  mixture.  The  authors 
realized  that  the  use  of  the  Fanning  expression  was  open  to  criticism  on 
the  ground  that  slip  was  neglected  and  that  the  fluids  lacked  homogeneity. 
In  the  region  of  high  gas-oil  ratios  studied  however,  these  assumptions 
were  not  as  serious  as  they  would  be  at  lower  rates. 

Rowels  (13)  studies  gas  lift  in  2 8  producing  wells  with  tubing  sizes 
varying  from  2-l/2  inches  to  6  inches.  Again,  data  were  plotted  as  a 
friction  factor  versus  a  "turbulence  factor",  equivalent  to  a  Reynolds 
number  using  average  input  mixture  properties. 

Martinelli  and  co-workers  (ll)  at  the  University  of  California 

developed  the  first  satisfactory  method  of  predicting  the  pressure  drop 

accompanying  the  simultaneous  flow  of  two  phase  mixtures  in  horizontal 

pipes.  Although  their  correlation  was  really  developed  for  horizontal 

flow,  the  method  of  approach  to  the  problem  is  worthy  of  mention.  These 

workers  developed  their  correlation  by  considering  four  flow  mechanisms 

(liquid  turbulent  -  gas  turbulent,  liquid  viscous  -  gas  turbulent,  liquid 

viscous  -  gas  viscous  and  gas  viscous  -  liquid  turbulent)  and  by  relating 

the  two  phase  pressure  drop  to  that  which  would  be  encountered  for  one  of 

the  phases  flowing  alone.  The  results  of  their  correlation  are  best 

expressed  mathematically: 

/APn  =  02  ,APx 

'AI/TP  Al/g 

where  /APs  =  the  pressure  drop  accompanying  two  phase  flow. 

W/TP 


'AP\  =  the  pressure  drop  with  the  gas  flowing  alone. 
AL  g 


- 
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0  =  a  function  of  X 

X  =  a  dimensionless  function  of  the  fluid  properties ,  tube 

diameter  and  hydraulic  radius,  different  for  each  flow  mechanism. 

In  a  paper  published  by  Lockhart  and  Martinelli  (10)  the  results  of 
all  preceding  investigations  were  recalculated  using  a  new  parameter  X  to 
replace  X  .  The  parameter  X  is  similar  to  %  but  it  is  more  general  and 
applies  to  all  flow  mechanisms.  The  parameter  X  is  defined  as  the  square 
root  of  the  ratio'  of  the  pressure  drop  in  the  pipe  if  the  liquid  flowed 
alone  to  the  pressure  drop'  if  the  gas  flowed  alone. 

In  1948  the  study  of  vertical  two  phase  flow  of  air  and  water  was 
started  at  the  University  of  Alberta  by  Radford  (15).  Radford  observed  the 
effect  of  water  and  air  flow  rates  on  the  flow  pattern  and  suggested  a  means 
of  calculating  the  pressure  drop  accompanying  two  phase  flow.  The  total 
pressure  drop  was  separated  into  the  hydrostatic  head  component  and  the 
irreversibility  component. 

Dunn  (4),  using  the  method  outlined  by  Radford,  correlated  the  pressure 
drop  in  terms  of  the  superficial  gas  phase  friction  factor  and  superficial 
gas  phase  Reynolds  number.  Tests  were  conducted  in  1.025  and  0. 63  inch 
diameter  tubes  at  an  average  pressure  of  l8  and  36  psia. 

Govier,  Radford  and  Dunn  (8),  applied  the  thermodynamic  analysis  to 
the  upward  vertical  flow  of  air -water  mixtures  in  a  1.025  inch  tube.  The 
pressure  drop  was  correlated  by  plotting  the  superficial  friction  factor, 
based  on  the  water  phase,  against  the  superficial  water  velocity.  The 
final  equations  developed  by  Govier,  Radford  and  Dunn  appear  below. 

Based  on  the  liquid  phase  -v  AP  1  +  R.m  1  f£F) 

L  AX  =  1  p  *  1  +  R  AX  L 

j-  .  Y 


(2) 


7 


(3) 


Alternately  for  the  gas  phase 


1  +  1/R^ 


1 


W  gcD 


2 

/AF\  = 


(5) 


where : 


P  =  absolute  pressure,  lbs.  per  sq.  ft.  absolute. 

AX=  tube  length,  feet. 

L  =  liquid  rate,  lbs-,  per  hr. 

G  =  gas  rate,  lbs.  per  hr. 

Vq=  specific  volume  of  gas,  cu.  ft.  per  lb. 

VL=  specific  volume  of  liquid,  cu.  ft.  per  lb. 

Em=  gas -liquid  mass  ratio,  G/L.  lbs.  per  lb. 

Rv=  gas-liquid  volume  ratio,  Gvq  ,  cu.  ft.  per  cu.  ft. 

Lvl 

Yjj=  superficial  liquid  velocity,  based  on  the  tube  cross 
section,  ft.  per  sec. 

Vq=  superficial  gas  velocity,  based  on  the  tube  cross 
section,  ft.  per  sec. 

D  =■  diameter,  feet. 

gc?=  dimensional  conversion  factor,  lb^ft 

lbRsec2 

f superficial  two  phase  friction  factor,  based  on  the  liquid 
phase ,  dimensionle s s . 


1q=  superficial  two  phase  friction  factor,  based  on  the  gas 
phase,  dimensionless. 


o 


Equations  (2)  and  (3)  were  employed  in  the  calculations  of  this 
work  hy  replacing  the  terms  applying  to  "liquid"  and  "gas"  with 
"reference  phase"  and  "second  phase"  respectively.  In  the  nomenclature 
used  here  equations  (2)  and  (3)  "become: 


ZiP  =  1  +  Rm  +  __1_  ,AS\ 
KSX  1  +  Kv  x  +.  Rv  '■AX''R 

.AF.  „  £4v|  ........... 

^/r  gcp 


(6) 


(7) 


The  "reference"  phase  is  defined  as  the  phase  having  the  lower  kinematic 
viscosity  -  water  in  all  cases  here. 

Short  (l6)  determined  the  effect  of  tubing  diameter  on  the  pressure 
drop  associated  with  the  upwards  vertical  flow  of  air-water  mixtures  in 
tube  diameters  of  O.63O,  1.025,  1.50,  and  2.50  inches.  The  mid-point 
pressure  was  kept  constant  at  3 6  psia.  Water  rates  varied  from  O.O695  to 
7.35  feet  per  second  and  air  rates  varied  from  0.OOO5H  to  O.523  cubic 
feet  per  second  at  test  section  pressure  and  temperature.  The  pressure 

1 

drop  was  correlated  by  plotting  the  superficial  friction  factor,  f p  , 
against  the  reduced  Reynolds  number  DVp  for  parameter  values  of  1)3  * 

Poettman  and  Carpenter  (l4)  developed  an  equation  for  the  flow  of 
gas-oil  mixtures  in  gas-lift  and  flowing  wells  using  the  mechanical  energy 
balance  and  Fanning  equations.  Density  of  the  mixture  was  calculated  from 
the  gas-oil  discharge  ratio.  Wo  attempt  was  made  to  calculate  the  viscosity 
as  it  was  assumed  to  have  a  negligible  effect  in  the  region  of  high 
turbulence  studied. 

In  1953  Galegar  (5)  of  the  University  of  Oklahoma  made  pressure  drop 
measurements  in  72  feet  of  2  inch  test  section  for  superficial  water 
velocities  from  0.1  to  1.6  ft.  per  sec.  and  superficial  air  velocities  from 
I.87  to  79.7  ft.  per  sec.  (measured  at  o0°F  and  14.7  psia).  Pressure  drops 


' 

’ 


■ 


9 


were  also  measured  in  19  feet  of  l/2  inch  pipe  for  superficial  velocities 

ranging  from  0.0125  to  O.96  ft.  per  sec.  for  the  water  and  from  I.85  to 

o 

28.8  ft.  per  sec.  for  the  air  (measured  at  60  F  and  14.7  psia).  The 
best  correlation  of  pressure  drop  data  was  found  to  be  a  plot  of  a 
modified  friction  factor,  based  on  the  input  gas  phase,  against  a  dim¬ 
ensionless  flow  ratio. 

Stovall  (18)  conducted  a  series  of  tests  on  the  vertical  flow  of 
kerosene-air  mixtures.  Equipment  constructed  by  Galegar  was  used  in 
the  experiments.  In  the  two  inch  test  section  the  maximum  superficial 
velocity  was  5 *6  ft.  per  sec.  for  the  kerosene  and  the  maximum  super¬ 
ficial  velocity  was  ^9  ft.  per  sec.  for  the  air,  (measured  at  standard 
conditions ) 

In  the  l/2  inch  test  section  the  maximum  superficial  kerosene 

velocity  was  1.23  ft.  per  sec.  and  the  maximum  superficial  air  velo- 

o 

city  was  37*5  ft.  per  sec.  (measured  at  60  F  and  14*7  psia).  Stovall 
correlated  his  data  using  Galegar* s  friction  factor  and  dimensionless 
flow  factor. 

A  paper  by  Galegar,  Stovall  and  Huntington  (6)  summarized  the 
data  obtained  by  Galegar  and  Stovall  by  presenting  a  single  plot  of 
the  friction  factor  versus  flow  ratio  for  both  the  kerosene-air  and 
air-water  systems. 

Equations  (6)  and  (7)  were  used  to  recalculate  the  basic  experi¬ 
mental  data  of  Galegar  (5)  and  Stovall  (18)  on  the  same  basis  as  the 
new  data  presented  in  Appendix  B-l.  Selected  values  are  reported  in 
Figures  30(a)  and  30(b). 
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In  1955  >  Calvert  and  Williams  (2)  analyzed  the  special  case  of 
upward,  vertical,  annular  flow  of  air-water  mixtures.  The  total 
pressure  drop  was  deined  in  terms  of  the  sum  of  two  factors,  the  pro¬ 
file  drag  and  the  skin  friction.  A  method  of  calculating  these  factors 
from  experimental  data  was  outlined,  and  a  comparison  between  cal¬ 
culated  and  experimental  results  was  presented. 

B.  FLOW  PATTERNS 

Observations  and  descriptions  of  flow  patterns  were  reported  in 
1936  and  1937  fry  Gosline  (7)  and  Cromer  and  Huntington  (3)  respect¬ 
ively.  These  preliminary  investigations  were  concerned  with  qualita¬ 
tive  considerations. 

Little  progress  was  made  in  relating  the  type  of  flow  pattern 
to  the  pressure  drop  until  Radford  (15)  described  the  flow  patterns 
with  respect  to  the  basic  pressure  drop  versus  air -water  volume 
ratio  curve. 

Govier,  Radford  and  Dunn  (8)  defined  four  pressure  drop  regimes 
in  terms  of  the  inflection  points  in  the  pressure  drop  curve.  These 
regimes  can  be  defined  precisely  whereas  flow  pattern  boundaries  are 
less  distinct.  Flow  patterns  were  described  (with  increasing  air 
rate  and  constant  water  rate)  as  bubble,  slug,  froth,  ripple,  film 
and  mist  flow. 

Calvert  and  Williams  (2)  noted  flow  patterns  similar  to  those 
observed  by  Govier,  Radford  and  Dunn  (8).  They  also  described  the 
flow  patterns  in  terms  of  the  pressure  drop  curve,  but  with  slightly 
different  flow  pattern  terminology. 
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C.  HOLDUP 

Moore  and  Wilde  (12)  and  Gosline  (7)  -were  the  first  invest¬ 
igators  to  recognize  the  importance  of  slippage.  Later,  Stein 
et  al  (17)  used  the  experimental  technique  employed  by  Gosline  (7) 
to  determine  the  gas  slippage.  A  semi-static  method  was  used  in 
which  gas  was  bubbled  through  a  column  of  oil.  The  average  gas 
velocity  or  slip  was  calculated  from  the  rise  of  liquid  in  the 
column  and  other  measurable  quantities. 

Dunn  (4),  Govier,  Radford  and  Dunn  (8)  and  Short  (l6)  invest¬ 
igated  slippage  in  the  upwards  vertical  flow  of  air-water  mixtures. 
Experiments  were  conducted  in  pipe  diameters  of  O.63  and  1.025  inches 
at  average  gas  pressures  of  18  and  36  psia.  The  term  holdup  ratio, 
used  here,  was  defined  by  Govier,  Radford  and  Dunn  (8)  as  the  dis¬ 
charge  air-water  ratio  divided  by  the  test  section  air-water  ratio. 
Slippage  is  related  to  the  holdup  by  the  following  equation  also 
presented  by  Govier,  Radford  and  Dunn  (8) : 


(8) 


» 


Where  =  average  velocity  of  the  gas,  ft/sec 


! 


=  avera£e  velocity  of  the  liquid,  ft/sec 


t  ! 


-  VL  =  slip  velocity,  ft/sec- 


In  terms  of  the  nomenclature  used  here: 
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III  Experimental  Equipment. 

A.  Air-Water  System. 

The  equipment  constructed  by  Radford  (15)  and  later  modified  by 
Dunn  (4)  and  Short  (l 6)  was  used  with  slight  alterations. 

The  test  section  comprised  32  feet  of  1.50  in.  seamless  copper 
tubing  equipped  with  full  flow  plug  valves.  Air  and  water 

were  admitted  at  the  bottom  of  the  test  section  through  a  standard 
l-l/2  in.  pipe  tee.  Water  was  supplied  from  the  main  at  45  psig  with 
provision  made  to  raise  this  pressure  to  220  psig  by  the  use  of  two 
5  H.P.  centrifugal  pumps  connected  in  series.  The  capacity  of  this 
system  was  70  I.G.P.M. 

A  suitable  quantity  of  air  was  available  from  the  laboratory  supply 
for  mid-point  pressure  of  18  psia  and  38  psia.  However,  at  pressures  of 
72  psia  and  110  psia  the  air  supply  had  to  be  supplemented  by  air 
compressors.  Two  rotary  compressors  with  a  capacity  of  125  S.C.F.M.  each 
and  a  maximum  pressure  of  200  psig  were  connected  in  parallel  directly 
into  the  air  line.  The  pressure  in  the  compressor  supply  tank  was  kept 
to  within  +  5  psig  by  a  pressure  control  valve  linked  directly  to  the 
engine  governor.  Air  from  both  the  laboratory  supply  and  compressor  was 
metered  with  an  orifice  plate  installation. 

The  air -water  mixture  was  discharged  from  the  top  of  the  test  section 
to  a  separator  where  a  constant  level  could  be  maintained  by  draining  the 
water.  Controlled  discharge  of  air  from  the  separator  kept  the  pressure  at 
the  mid  point  of  the  test  section  constant.  The  air  discharged  from  the 
separator  through  a  pneumatic  valve  operated  by  a  0  -  100  psig  Brown 
pneumatic  pressure  controller.  The  pneumatic  controller  was  activated  by 
the  mid-point  pressure  signal. 
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The  total  pressure  drop  was  measured  over  22.88  feet  of  the  test 
section  with  the  lower  tap  approximately  8  feet  above  the  mixing  tee. 

Two  60  inch  manometers  containing  water  and  mercury  respectively  were 
used  to  measure  the  pressure  drop.  A  single  phase,  air  from  small  plastic 
separators  located  at  the  two  pressure  stations,  was  employed  as  the  lead 
line  fluid.  Provision  was  made  to  measure  the  air  and  water  temperatures. 

B.  Oil-Water  System. 

Extensive  modifications  were  made  to  the  air-water  system  equipment 
as  shown  in  the  schematic  diagram.  Figure  1. 

The  oil-water  test  section  consisted  of  32  feet  of  1.025  inch  I.D. 
cellulose  acetate  butyrate  tubing.  It  was  equipped  with  Nordstrom  full 
flow  plug  valves  -  similar  to  the  air -water  system. 

Water  and  oil  were  pumped  simultaneously  through  the  test  section  to 
a  gravity  separator.  The  oil  was  decanted  to  a  45  gallon  storage  container 
while  the  water  was  drained  from  the  bottom  of  the  separator  to  a  similar 
container. 

Oil  from  the  storage  container  was  pumped  through  the  system  by  a 
"builtogether15  5  H.P.  centrifugal  pump  capable  of  delivering  70  gallons 
per  minute  with  a  shut  off  head  of  80  psig.  The  oil  was  metered  through 
one  of  three  rotameters  covering  a  range  of  flow  from  0  to  25  I.G.P.M. 

A  similar  pump  was  used  for  the  water.  Three  rotameters  were  used  to 
measure  the  water  rates  from  0-25  I.G.P.M. 

The  temperatures  of  the  oil  and  water  streams  were  measured  with  dial 
type  bimetallic  thermometers  mounted  in  the  discharge  lines  from  the 
rotameters. 

The  total  pressure  drop  was  measured  over  22.75  feet  of  the  test 
section  with  water  over  mercury  and  water  over  carbon  tetrachloride 
manometers  connected  in  parallel.  Single  phase  lead  line  fluid  was 
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PIG.  1  Schematic  Plow  Diagram  of  Equipment 


15 


obtained  from  plastic  oil-water  separators  located  at  the  two 
pressure  taps. 
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IV  Experimental  Procedure . 

A.  Air- Water  System. 

The  tests  were  conducted  using  the  same  procedure  as  described 
by  Short  ( l6 ) . 

Mid  point  pressures  of  18,  72  and  110  psia  were  investigated  at  a 
number  of  constant  water  rates  in  the  1.5  inch  copper  tube.  With  the 
water  rate  set  at  a  constant  value ,  the  air  rate  was  varied  from  the 
lowest  to  the  highest  value  at  which  the  mid  point  pressure  could  be 
controlled. 

A  typical  test  consisted  of  setting  a  fixed  water  rate  and  an 
arbitrary  air  rate.  When  steady  state  had  been  reached  the  air  and 
water  rates,  mid-point  pressure  and  temperatures  were  recorded.  The 
full  flow  valves  were  then  closed  simultaneously  and  the  volume  of  water 
in  the  test  section  was  measured. 

B.  Oil-Water  System. 

The  series  of  tests  were  carried  out  at  8  different  constant  water 
rates  in  the  1.025  inch  plastic  tube.  The  oil  rate  was  varied  from  zero 
to  the  maximum  for  each  water  rate. 

A  typical  test  consisted  of  setting  a  fixed  water  rate  on  the  rotameter 
and  an  arbitrary  oil  rate.  The  oil  water  interface  was  adjusted  to  a 
constant  level  in  the  separator.  When  steady  state  was  reached  the 
pressure  drop,  oil  and  water  rates,  oil  and  water  temperatures  and  room 
temperature  were  recorded.  As  the  second  phase  fluid,  oil,  was 
incompressible  it  was  not  necessary  to  record  the  mid-point  pressure. 

The  holdup  was  then  obtained  by  closing  the  full  flow  plug  valves, 
allowing  the  oil  and  water  to  separate,  and  reading  the  volume  of  oil 


on  the  calibrated  tube. 
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At  low  flow  rates  visual  observations  of  the  flow  patterns  were 
recorded.  At  higher  rates  Superficial  water  velocities  in  excess  of 
1  ft .  per  sec . )  photographs  were  obtained  to  show  the  type  of  flow 


pattern. 


V  Experimental  Data 

The  basic  experimental  data  are  presented  in  tabular  form  in 
appendices  A  and  C. 

Table  C-l  of  Appendix  C  gives  the  calibration  data  for  the  water 
and  oil  rotameters.  The  data  are  reported  as  mass  rates  of  flow  in 
lbs.  per  sec.  at  even  increments  of  the  rotameter  scale.  The  data 
are  also  represented  graphically  in  Figures  Cl-G 6  (Appendix  C)  as 
mass  rates  of  flow  divided  by  rotameter .setting  versus  rotameter 
setting. 

Figures  2  to  17  are  constructed  from  data  presented  in  columns 
5,  7  and  8  of  Table  A-2  (Appendix  A).  These  figures  show  the  total 
unit  pressure  drop  and  holdup  ratio  as  functions  of  a  common  abscissa,, 
the  air-water  volume  ratio.  The  pressure  drop  and  holdup  ratio  curves 
are  similar  in  shape  to  those  obtained  by  Short  (l6)  and  Dunn  (4). 

The  corresponding  pressure  drop  curves  for  the  oil-water  system 
(Figures  18  -  24)  are  obtained  from  data  tabulated  in  Table  A- 3  of 
Appendix  A.  The  oil-water  pressure  drop  curves  resemble  the  air-water 
pressure  drop  curves  at  superficial  reference  phase  velocities  above 
one  foot  per  second.  However,  below  one  foot  per  second  there  is  no 
evidence  of  a  maximum  or  second  minimum  in  the  curves. 

The  pressure  drop  curves  (Figures  2  -  24)  are  divided  into  two  or 
more  regimes  according  to  the  method  of  Govier,  Radford  and  Dunn  (8). 
Points  on  the  pressure  drop  curves  corresponding  to  the  first  minimum, 
maximum  and  second  minimum  locate  the  termination  of  regimes  I,  II  and 
III  respectively. 
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In  Figure  25  the  loci  of  the  first  minimum  of  the  pressure  drop 
curves  are  plotted  as  a  function  of  the  superficial  reference  phase 
velocity  for  each  of  the  second  phase  densities.  These  loci  are 
parallel  straight  lines  on  logarithmic  paper  for  all  air  densities 
over  the  range  of  superficial  reference  phase  velocities  studied. 

The  locus  of  the  first  minimum  for  the  oil-water  system  is  similar  to 
that  of  the  air -water  system  for  superficial  reference  phase  velocities 
up  to  one  foot  per  second.  Above  one  foot  per  second  the  locus  curves 
towards  lower  values  of  R  than  would  be  obtained  from  a  continuation 
of  the  linear  relationship. 

In  Figure  2 6  the  loci  of  the  maximum  and  second  minimum  of  the 
pressure  drop  curves  are  plotted  against  the  superficial  reference 
phase  velocity.  With  air  as  the  second  phase  the  loci  of  the  maximum 
of  the  pressure  drop  curves  are  parallel  straight  lines  on  logarithm 
paper  up  to  a  superficial  reference  phase  velocity  of  1.5  feet  per 
second.  Above  this  value  the  maximum  loci  curve  towards  the  second 
minimum  loci  until,  at  a  superficial  reference  phase  velocity  of  5*5 
feet  per  second,  the  maximum  and  second  minimum  become  coincident. 

It  was  not  possible  to  determine  this  point  of  convergence  experi¬ 
mentally  due  to  an  insufficient  supply  of  air. 

When  plotted  against  the  superficial  reference  phase  velocity 
the  second  minimum  loci  for  the  air  are  also  parallel  straight  lines 
terminating  at  the  point  of  convergence  with  the  maximum.  With  oil 
as  the  second  phase,  however,  both  the  maximum  and  second  minimum 
become  coincident  with  the  first  minimum  at  superficial  reference 
phase  velocities  below  one  foot  per  second.  Observations  in  the  lab- 
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oratory  and  photographs  help  to  clarify  this  statement. 

At  low  superficial  reference  phase  velocities  the  transition  from 
bubble  and  slug  flow  of  regime  I  to  mist  flow  of  regime  IV  occurs  in 
a  narrow  region  on  the  pressure  drop  curve.  This  is  illustrated  with 
reference  to  the  photographs  of  Figure  27  and  the  corresponding  pressure 
drop  curve.  Figure  18.  Part  A  of  Figure  27  represents  bubble  flow  at 
a  superficial  reference  phase  velocity  of  0.1  feet  per  second  and  a 
volume  ratio  of  approximately  unity.  Observations  in  the  laboratory 
indicate  bubbles  of  oil  in  a  continuous  phase  of  water.  As  the  oil 
rate  is  increased  the  oil  bubbles  cluster  together  to  form  slugs 
(Fig.  27-B).  The  leading  edge  of  the  slugs  is  not  clearly  defined  as 
.in  the  case  of  air-water  flow.  The  oil  slugs  tend  to  be  shorter  than 
the  air  slugs  and  are  completely  surrounded  with  bubbles  of  oil  which 
makes  them  harder  to  distinguish.  As  the  oil  rate  is  further  increased 
the  pattern  changes  to  a  condition  corresponding  to  mist  flow  in  the 
air  water  system.  Water  droplets  are  entrained  in  a  continuous  oil 
phase.  This  condition  is  illustrated  in  Figure  27-C  in  which  droplets 
of  water  can  be  seen  in  the  oil.  Observations  in  the  laboratory  sub¬ 
stantiate  the  fact  that  oil  is  the  continuous  phase. 

Figure  28  illustrates  the  flow  pattern  at  various  volume  ratios 
at  a  superficial  reference  phase  velocity  of  1  ft.  per  sec.  (See  Fig¬ 
ure  20).  Again  bubble  flow  (oil  in  water)  is  evident  at  a  low  |K 
(Figure  28  A) .  As  the  oil  rate  is  increased  a  pattern  between  slug 
flow  and  froth  flow  is  apparent  from  Figure  28B.  Finally,  in  Figure 

28  C  at  an  R  of  5*24  a  condition  of  mist  flow  is  reached.  It  is 
V 

not  too  clear  from  the  photograph  which  is  the  continuous  phase. 
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However,  it  was  observed,  by  stopping  the  flow,  that  small  bubbles  of 
water  descended  through  a  continuous  oil  phase  as  opposed  to  oil  rising 
through  the  water  in  the  tests  represented  by  Figures  28  A  and  28  B. 

The  photographs  of  Figure  29  are  taken  at  a  superficial  reference 
phase  velocity  of  2.6  feet  per  second  where  the  pressure  drop  curve 
(Figure  22)  exhibits  a  maximum  and  second  minimum.  Figure  29  A  was 
taken  at  an  R^  value  of  0.48  where  the  beginning  of  slug  flow  is 
evident.  In  Figure  29  B  the  flow  pattern  is  definitely  froth  flow. 

This  flow  pattern,  not  distinguishable  at  lower  flow  rates,  is  char¬ 
acterized  by  a  thorough  mixing  of  the  phases.  A  longer  time  is  required 
for  complete  separation  of  the  phases  In  such  a  flow  pattern  than  at 
.higher  R^  values  where  the  oil  is  a  continuous  phase.  Photographs  In 
the  mist  flow  region  at  a  superficial  reference  phase  velocities  of 
2.6  ft  per  sec.  were  blurred  and  are  not  included. 

The  holdup  ratio  curves  for  the  air-water  system  are  similar  in 
shape  to  those  obtained  by  Short  (l6)  and  Dunn  (4).  However,  a  marked 
decrease  in  holdup  ratio  is  evident  in  the  oil-water  system;  the  ratio 
being  close  to  unity  for  all  superficial  reference  phase  velocities. 

A  scattering  of  the  data  is  noticeable  at  low  volume  ratios  for  all 
superficial  reference  phase  velocities.  A  possible  explanation  lies 
in  the  fact  that  there  was  a  large  percentage  error  in  reading  the 
holdup  ratio  with  the  small  amount  of  oil  in  the  tube  at  low  oil-water 


ratios. 
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VI  Interpretation  of  Results 

The  following  is  a  discussion  of  the  steps  taken  to  process  the 

experimental  data  and  to  determine  the  effect  of  second  phase  density 

on  the  pressure  drop,  flow  patterns  and  holdup. 

(a)  The  Effect  of  Second  Phase  Density  on  Pressure  Drop. 

Equations  (6)  and  (7)  were  used  to  calculate  the  superficial 

friction  factor,  f^,  based  on  the  superficial  reference  phase  velocity. 

As  explained  further  in  Appendix  D  the  friction  factor  was  plotted 

as  a  function  of  the  volume  ratio  to  form  a  ba,sic  "working  plot" 

(Figure  D-l).  From  a  complete  set  of  these  "working  plots"  -  one 

for  each  of  the  densities  and  diameters  studied  -  a  series  of  cross 

plots  were  prepared  leading  to  Figure  D-8  (Appendix  D).  This  is  a 

nlot  of  f’(DV  (vo  -  is )~^*33  versus  DJ’^  V  and  indicates  that 

R  R  s  s 

fp  (DVj^5  (l/Q  -V  )  depends  upon  the  1.1  power  of  D"*^Vs* 

The  difference  in  kinematic  viscosity  of  the  phases  was  chosen, 
after  many  trials,  as  the  best  way  of  expressing  the  effect  of  den¬ 
sity  and  (for  the  values  covered)  of  viscosity. 

In  Figures  30-a  and  30-b  the  relationship  between  f 1  and  DV 

0  ^  R  ** 

is  represented  for  a  range  of  D0*°  V^1*1  {l/Q  ~V^)  parameters 

from  0.0001  to  0.1.  A  range  of  f 1  values  from  0.00^-  to  1000  is 

R 

7 

covered  in  Figure  30-a  and  continued  to  an  f'  of  2  x  10  in  Figure 

X\ 

30-b.  In  order  to  verify  the  correlation  and  test  for  errors  in- 

O  Oj- 

troduced  by  smoothing,  a  complete  set  of  values  of  f,  DVp  and  v'  ^ 

„  1.1,,,  \°-31  R 

vs  \  were  calculated  from  the  basic  data.  Da/fca  reported 

here  (Appendix  B)  were  supplemented  with  data  collected  by  Short  (l6) 


V  , 


If  ' 

. 

. 
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and  Dunn  (k) .  Additional  data  were  selected  from  the  kerosene-air 

system  investigated  by  Stovall  (l8),  the  air-water  system  studied 

by  Galegar  (5)  and  the  air-lift  pump  tests  conducted  by  Davis  and 

Weidner  (l)  at  the  University  of  Wisconsin. 

Where  the  f '  versus  R  relationship  was  readily  available  a 
R  v 

number  of  points,  necessary  to  define  the  curve,  were  chosen  and 
the  value  of  D^*u^  (^s  "*z'p)  was  calculated  for  each, 

point.  In  the  case  of  the • data  obtained  by  Dunn,  Galegar,  Stovall 
and  Davis  and  We idner,  a  number  of  individual  experimental  points 
were  chosen  at  random  and  calculated  similarly.  These  points  appear 
on  the  final  correlation  (Figures  30~a  and  3^-b)  with  numerical 
superscripts  equal  to  I?*®  Vg1*^  V  g-  ^p)°* 33  and  alphabetical  sub¬ 
scripts  referring  to  the  tube  diameter. 

The  Hagen- Foiseuille  and  Karman-Nikur ad se  lines  for  single 
phase  fluids  are  shown  on  Figure  30_a  for  reference  purposes. 

The  superficial  Reynolds  number  based  on  the  reference  phase  may 
be  obtained  by  multiplying  the  abscissa,  DV^,  by  Pp/|j.-p>,  the  numeri- 

5  0 

cal  value  of  which  is  1.02  x  10  at  75  F.  The  ordinate  f^  becomes 

f  when  the  system  is  single  phase.  The  two  phase  friction  factor, 

f’,  is  different  from  a  single  phase  f  for  two  reasons.  Firstly, 

R 

the  two  phase  friction  factor  is  based  on  a- superficial  reference 

phase  velocity  which  increases  the  f 1  by  a  function  of  the  volume 

R 

ratio.  Secondly  the  two  phase  friction  factor  accounts  for  all 
irreversibilities  of  the  two  phase  system.  This  effect  is  super¬ 
imposed  on  the  increase  due  to  the  volume  ratio  with  the  net  effect 
that  the  two  phase  friction  factors  always  lie  above  the  single  phase 


line. 


0.002  0005  0.01  0.02  0.05 


As  the  reference  phase  velocity  is  increased  at  a  fixed  value  of 
DJ*^V 


J (vo  a  point  is  reached  where  the  linear  logarithmic 


R- 

relationships  between  f ’  and  V  does  not  hold.  The  line  for  a  given 
parameter  value  curves  and  eventually  blends  into  the  single  phase 
Karman-Nikuradse  line  at  a  higher  DV^  as  shown  in  Figure  30-a  for 
parameter  values  of  0.00001  and  0.00003*  At  lower  superficial  reference 
phase  velocities  the  logarithmic  relationship  between  f^  and  DV^  remains 
linear.  This  is  verified  to  values  of  DV^,  as  low  as  0.002  (f^  values 

y 

as  high  as  10  ) . 

Figures  30-a  and  30 -Id  illustrate  the  effect  of  the  second  phase 
density  on  the  superficial  friction  factor  and  therefore  directly  on  the 
pressure  drop  due  to  irreversibilities.  Reference  is  made  to  basic 
equations  (6)  and  (7)  for  the  effect  of  second  phase  density  on  the  total 
unit  pressure  drop  (irreversibilities  +  hydrostatic  head). 

Substituting  equation  (7)  into  equation  (6)  gives 


VR  ^  +  ^m  + 


AX 


with 


1  +  R-< 


R 


1  +  *, 

constant 

constant 

constant 


(ISA ) 

gcD 


(10) 


L  =  constant 
V  =  constant 


equation  (10)  becomes 

vR  AP  =  1  +  psj[D2/4  + 

**  1  +  ^PrA 


1  +  rcD^P^/  4L  (■ 


2  f  R  VR 


ec- 


(11) 


Increasing  the  second  phase  density,  ps,  increases  the  hydrostatic  head 

component  and  thus  the  total  unit  pressure  drop.  The  irreversibility 

component  on  the  other  hand  is  decreased  with  increasing  ps  (Figures  30-a 

and  30-b).  The  net  effect  on  the  total  pressure  drop  depends  on  the 
relative  quantity  of  second  phase  to  "reference  phase. 
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(B)  The  Effect  of  Second  Phase  Density  on  the  Flow  Patterns. 

The  effect  of  second  phase  density  on  the  flow  patterns  is  reflected 
in  the  loci  of  the  first  minimum,  second  minimum  and  maximum  of  the 
pressure  drop  curves.  In  general,  an  increase  in  the  second  phase 
density  shifts  the  loci  of  the,  first  minimum,  maximum  and  second 
minimum  to  a  lower  value  of  the  volume  ratio  at  a  given  superficial 
reference  phase  velocity  as  shown  in  Figures  25  and  26. 

The  first  minimum  in  the  pressure  drop  curve  is  correlated  in 

Figure  31  for  all  second  phase  densities  as  a  plot  of  superficial 

-0.5 

reference  phase  velocity  versus  *  Again,  kinematic  viscosity 

is  chosen  as  the  best  way  of  expressing  density  and  (for  the  values 
covered)  viscosity.  The  correlation  is  obtained  from  Figure  25  by 
determining  the  relationship  between  the  volume  ratio,  R^,  and  the 
kinematic  viscosity,!/,;  at  constant  superficial  reference  phase  velocity. 
The  volume  ratio,  R^,  is  found  to  vary  with  V® *  ^  at  constant  V  .  The 

first  minimum  of  the  pressure  drop  curve  is  not  a  function  of  the 

-0.5 

diameter  as  shown  by  Short  (l6).  The  VR  versus  ^VQ  *  relationship 

(Figure  31 )  is  linear  for  all  second  phase  densities  up  to  a  superficial 

reference  phase  velocity  of  one  foot  per  second.  Above  this  velocity 

the  non-linearity  of  the  basic  data  (Figure  25)  for  the  oil  is 

evident.  This  non-linearity  may  be  noticeable  in  the  air-water  system 

at  higher  superficial  reference  pha.se  velocities.  If,  however,  it  is 

an  absolute  viscosity  dependent  function,  studies  now  underway  at  this 

University  will  show  this  to  be  so. 

In  Figure  32  the  second  minimum  of  the  pressure  drop  curve  is 

correlated  by  plotting  the  superficial  reference  phase  velocity,  VR, 

against  the  quantity  R  p  ^ ^ .  The  relationship  between  R  and  p 

vs  vs 

is  found  from  Figure  26  by  plotting  Rv  against  pg  at  constant  VR.  Again, 
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the  effect  of  diameter  is  obtained  from  Short  (l6)  who  determined 
the  quantity  0  by  trial.  The  locus  of  the  second  minimum  of  the 
pressure  drop  curve  for  the  oil-water  system  does  not  appear  on  the 
correlation  (Figure  32)  because  there  is  insufficient  overlapping  of 
the  oil-water  and  air-water  second  minimum  loci  to  justify  their 
correlation.  Below  a  superficial  reference  phase  velocity  of  one 
foot  per  second  the  maximum  of  the  pressure  drop  curve  for  the  oil 
disappears  and  the  first  minimum  and  second  minimum  are  coincident 
(Figure  26).  This  suggests  that  the  locus  of  the  first  minimum  of 
the  pressure  drop  cruve  for  the  oil  (Figure  25)  will  also  fit  the 
second  minimum  correlation  below  a  VR  of  1  ft.  per  sec.  A  calculation 
at  several  points  reveals  that  this  is,  in  fact,  the  case. 

Figure  33  represents  the  locus  of  the  maximum  of  the  pressure  drop 

curve  as  a  function  of  the  superficial  reference  phase  velocity  and  the 

quantity  R  l/2.  Agapn  the  quantity  was  obtained  from 

V  s 

Short  (l6).  The  Rv  -  Ps  relationship  is  determined  from  Figure  26  by 
plotting  Ry.  against  pg  at  a  constant  The  maximum  converges  with 

the  second  minimum  as  shown  by  Figure  26.  This  convergence  is  a  function 
of  dipmeter  as  shown  by  Short  (l6).  Data  of  Short's  for  a  second  phase 
density  of  0.184  lbs.  per  cu.  ft.  at  diameters  of  O.63O  ins.  and  1.025  ins. 
are  plotted  on  Figure  33.?  although  no  data  exists  for  other  second  phase 
densities  at  these  diameters.  The  convergence  points  are  indicated  by 
dashed  lines  on  Figures  32  and  33-  As  the  maximum  of  the  pressure  drop 
locus  for  the  oil  does  not  exist  at  superficial  reference  phase  velocities 
below  1  foot  per  second,  it  cannot  be  correlated  in  this  region.  There¬ 
fore  the  correlation  of  the  loci  of  the  maximum  of  the  pressure  drop 
curves  is  limited  to  second  phase  densities  up  to  0.552  lbs.  per  cu.  ft. 
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(c)  The  Effect  of  Second  Phase  Density  on  Holdup. 

Figure  39,  a  plot  of  H^V^  °*00^2  Pr  ps ^ versus  pR  -  at 

l/o 

constant  R^D  '  J ,  shows  the  effect  of  second  phase  density  on  holdup 
for  regimes  I,  II  and  III.  The  holdup  ratios  from  Tables  A-2  and  A- 3 
are  plotted  in  Figures  3^  to  38  as  holdup  ratio ,  HR,  versus  superficial 
reference  phase  velocity,  VR,  at  constant  B  (The  quantity  RVD‘L///^ 

is  obtained  from  work  presented  by  Short  (l6)  and  is  valid  in  pressure 
drop  regimes  I,  II  and  III.)  In  addition  data  of  Short’s  for  a  second 
phase  density  of  0.184  lbs.  per  cu.  ft.  is  plotted  on  Figures  34  to  38  • 
There  is  little  difference  in  the  holdup  ratio  when  air  is  the 
second  phase.  However,  with  oil  as  the  second  phase  the  holdup  ratio 
is  much  smaller,  being  near  unity  for  all  superficial  reference  phase 
velocities.  This  suggests  that  the  holdup  ratio  is  some  function  of  the 
difference  between  the  reference  phase  density  and  second  phase  density. 
A  best  fit  line  is  then  put  through  the  data  plotted  in  Figures  3^  to  38 
for  ( Pp  -  pg)  values  of  approximately  62  lbs.  per  cu.  ft.  for  the  air 
and  l4  lbs.  per  cu.  ft.  for  the  oil.  Since  HR  is  a  function  of  VR 
when  (p^  -  ps)  is  approximately  62  lbs.  per  cubic  ft.,  and  is  obviously 
not  a  function  of  YR  when  (pj>  -  pg)  is  zero,  the  way  in  which  the 
%  -  %  relationship  varies  with  density  is  determined. 

Stated  mathemat i c ally : 

m 

%  a  VR 

where  m  =  0 1 ( p^  -  ps ) 

The  power  m  is  determined  by  plotting  m  versus  (pj^  -  p  )  for  air,  oil 

and  water  density  differences.  A  linear  relationship  between  rn  and 

(pp  -  P  )  is  obtained,  m  being  equal  to  0.0022(p  -  p  ).  The  final 
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plot  (Figure  39)  is  obtained  by  plotting  Hp>V-^  ps)  versug 

(p^  -  ps)  for  values  of  Qf  1,  2,  5;  10  and  20  obtained  from 

Figures  3^-  to  3$* 

It  can  be  seen  from  Figure  39  that  an  increase  in  the  second  phase 
density  decreases  the  holdup  ratio.  Also,  as  the  diameter  decreases, 
the  holdup  ratio  decreases,  which  is  consistent  with  capillary  flow 


theory. 


: 


. 


wf 


68 


VII  Conclusions 

(l)  The  upwards  vertical  flow  of  two  phases,  a  second  phase  of 
air  or  oil  and  a  reference  phase  of  water,  has  been  investigated  for 
the  following  range  of  variables: 


Tube 

Second 

Second  Phase 

Second  Phase  Second  Reference 

Diameter 

Phase 

Absolute 

Kinematic  Phase  Phase 

ins. 

Density 
lbs/ ft0 - 

Viscosity* 
lbs/ft. sec. 
x  105 

Viscosity*  Rates  Rat^s 

ft^/sec.  ft^/sec.  ft. /sec. 

x  io5 

1.50 

0.092 

1.21 

13.1  0.021-0.057  0.000851-0.01065 

1.50 

O.362 

1.21  • 

3.36  0.024-0.403  0.000851-0.0901 

1.50 

0.552 

1.21 

2.20  0.0215-0.225  0.000851-0.0425 

1.025 

48.57 

62.4 

1.29  0.00026-0.0677  0.000573-0.0322 

(2)  The  friction  factor,  f.^  has  been  correlated  with  second  phase 
velocity,  reference  phase  velocity,  tube  diameter  and  the  difference 
in  kinematic  viscosity  of  the  phases. 

(3)  The  pressure  drop  accompanying  two  phase  flow  can  be  calculated 
for  any  second  phase  density  from  Figures  30-a  and  30-b  with  the  aid 

of  equations  (6)  and  (7) 

(4)  The  f^  correlation  shows  agreement  with  data  from  the  litera- 
ture  (l),  (5),  (18). 

(5)  The  flow  patterns  have  been  correlated  in  terms  of  the  loci 
of  the  inflection  points  in  the  basic  pressure  drop  curves. 

(6)  The  first  minimum  in  the  pressure  drop  curve  has  been  correl- 
ated  in  Figure  31  for  all  second  phase  densities. 

(7)  The  second  minimum  in  the  pressure  drop  curve  can  be  obtained 
from  Figure  3 2  for  second  phase  densities  from  0  to  0.552  lbs.  per 


cu.  ft. 


*  at  test  section  conditions. 
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(8)  In  Figure  33  the  maximum  of  the  pressure  drop  curve  has  been 
correlated  for  second  phase  densities  from  0  to  0.552  lbs.  per  cu.  ft. 

(9)  The  holdup  ratio  for  regimes  I.,  II  and  III  can  be  obtained 
from  Figure  39  "by  knowing  the  tube  diameter,  second  phase  velocity 
and  second  phase  density. 

(10)  The  limiting  condition  that  the  holdup  ratio  approach  unity 
for  a  single  phase  fluid  is  satisfied  in  Figure  39* 

(11)  When  studies  of  the  effect  of  second  phase  viscosity  are 
completed,  generalized  correlations  of  the  pressure  drop,  flow  patterns 
and  holdup  for  all  fluid  two  phase  systems  will  be  available. 
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P 


T 


V 

L 

V 
R 


V 


s 


Vf 

G 


NOMENCLATURE 

=  pipe  diameter,  ft. 

=  superficial  two  phase  friction  factor,  based  on  the  gas  phase 

=  superficial  two  phase  friction  factor,  based  on  the  liquid  phase 

=  superficial  reference  phase  friction  factor,  dimensionless 

=  superficial  second  phase  friction  factor,  dimensionless 

=  saturated  second  phase  rate,  lb/sec 

=  dimensional  conversion  factor,  lbm  ft 

Ibw  sec2 

=  holdup  ratio,  dimensionless 

=  reference  phase  rate,  lb/sec 

=  mass  ratio,  G  ,  lb.  air  Qr  lb  Qil 

L  '•lb.  water  1^' water 


=  volume  ratio  G 

L  R 

=  pressure,  lb/ft2  absolute 
=  temperature,  °F 

=  superficial  gas  velocity,  based  on  the  tube  cross  section,  ft/sec. 

=  superficial  liquid  velocity,  based  on  the  tube  cross  section, ft/sec. 
=  superficial  reference  phase  velocity,  based  on  the  cross  sectional 
area  of  the  tube,  ft/sec. 

=  superficial  second  phase  velocity,  based  on  the  cross  sectional 
area  of  the  tube,  ft. /sec. 

=  average  velocity  of  the  gas,  ft/sec. 

=  average  velocity  of  the  liquid,  ft/sec. 


.3 


air 


’ft ^  water 


or 


S  . 

ft  air 
ft3  water 


or 


ft3oil 
ft 3  water 
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V*  =  average  lineal  velocity  of  the  reference  phase,  ft. /sec. 
R 

V’  =  average  lineal  velocity  of  the  second  phase,  ft/sec. 
s 

v  =  specific  volume  of  the  reference  phase,  ft.^/lb. 

R 

v  =  specific  volume  of  the  second  phase,  ft.^/lh. 


AX 


tube  height,  ft, 


y\p 

=  the  pressure  drop  acoompanying  two  phase  flow. 


V\L  'TP 


the  pressure  drop  with  gas  flowing  alone. 


Greek  Letters: 


V  =  kinematic  viscosity  of  the  reference  phase,  ft.  /sec. 

R 

7/  =  kinematic  viscosity  of  the  reference  phase,  ft.  /sec. 

.  ^  s 

u  =  absolute  viscosity  of  the  reference  phase,  lb/ft.  sec, 
R 

jj.  .=  absolute  viscosity  of  the  second  phase,  lb/ft.  sec. 
s 

o 

p  =  density  of  the  reference  phase,  lb/ft. 
p  =  density  of  the  second  phase,  lb/ft^ 


s 

0  =  a  function  of  the  dimensionless  quantity 

*  =  indicator  of  a  functional  relationship 
JC  =  a  dimensionless  function  of  the  fluid  properties,  tube  diameter 
and  hydraulic  radius. 

Subscripts 
G  =  gas  phase 
L  =  liquid  phase 

R  =  reference  phase,  defined  as  the  phase  with  the  lower  kinematic 
viscotity 


s  =  second  phase 


1 


72 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 

(12) 

(13) 

(14) 
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APPENDIX  A 

1.  Supplementary  Data. 

2.  Experimental  Data  -  Air-Water  System. 

3*  Experimental  Data  -  Oil-Water  System. 


tllfe 
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APPENDIX  A 


Table  A-l  contains  data  on  the  geometry  of  the  oil-water  and 
air-water  systems.  In  addition,  the  physical  properties  of  the  oil 
are  listed. 

The  basic  experimental  data  for  the  air -water  system  are  presented 
in  columnar  form  in  Table  A-2. 


Column 

1 

2  ‘ 

3 

k 

5 

6 

7 


DATA 

Test  number. 

Water  rate  in  cubic  feet  per  second  calculated 
at  mid-point  temperature. 

Saturated  air  rate  calculated  at  mid -point  pressure 
and  temperature. 

Mid-point  temperature,  °F. 

Discharge  air-water  ratio,  cu.  ft.  of  saturated 
air  per  cu.  ft.  of  water,  obtained  by  dividing 
column  3  hy  column  2. 

Test  section  air-water  ratio,  obtained  by  measuring 
the  volume  of  water  and  calculating  the  air  volume 
by  difference. 

Holdup  ratio,  defined  as  the  discharge  air-water 
ratio  divided  by  the  test  section  air-water  ratio 
(column  5  divided  by  column  6). 

Unit  pressure  drop,  feet  of  water  per  foot  of  tube. 


8 


obtained  by  converting  basic  pressure  drop  data  to 
feet  of  water  and  dividing  by  the  test  section  length. 


' 


■ 


... 


■ 
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Table  A-3  contains  experimental  data  for  the  oil-water  system. 
The  columns  in  Table  A-3  are  equivalent  to  those  for  the  air-water 
system  (Table  A-2)  described  above  when  the  word  "air"  is  replaced 


with  "oil". 


V 


m 
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TABLE  A-l 
Supplementary  Data 


(1)  Air -Water  System 

Tube  diameter  =  1.50  ins. 

Test  section  length  =  22.88  feet. 

Volume  between  plug  valves  =  10,950  c.c. 

(2)  Oil-Water  System 

Tube  diameter  =  1.025  ins. 

Test  section  length  =  22.75  feet. 

Volume  between  plug  valves  =  5*200  c.c. 

Oil  density  =  48.57  lbs.  per  cu.  ft.  at  75°F. 

Oil  viscosity  =  1.20  centistokes  at  75°F 

=  1.29  x  10 sq.  ft.  per  sec.  at  75°F. 


EXPERIMENTAL  DATA 
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APPENDIX  B 


1.  Calculated  Data  -  Air-Water  System 


2. 


Calculated  Data  -  Oil-Water  System 
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APPENDIX  B 

Tables  B-l  and  B-2  contain  calculated  results  of  the  experimental 
data  for  the  air-water  and  oil-water  systems  respectively. 


Column 

1 

2  ' 


Data 


Test  Number 

Total  Unit  Pressure  drop,  v^>  AP,  feet  of  water  per 

AX 

foot  of  tube. 

Hydrostatic  head  component  of  the  total  pressure 
drop,  1  +  Rffi  ,  feet  of  water  per  foot  of  tube. 

1  +  Ry 

Irreversibility  component  of  the  total  pressure  drop, 

AF\  ,  feet  of  water  per  foot  of  tube  (column  2 


(%) 


1  +  Ry  AX'r 
minus  column  3)* 

Superficial  reference  phase  velocity,  feet  per  second, 
where  the  reference  phase  is  defined  as  the  phase 
having  the  lower  kinematic  viscosity  (water). 

Two  phase  friction  factor,  f-^,  based  on  the  super¬ 
ficial  reference  phase  velocity,  defined  by  the 


relation  (&F) 


=  2  f*  V 

XR  VR 

ScD 


Pressure  drop  regime 


TABLE  B-l 
CALCULATED  DATA 
AIR-WATER  SYSTEM 
TUBE  DIAMETER  1.50  INS 
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APPENDIX  C 
Rotameter  Calibrations 


The  six  rotameters  were  calibrated  by  a  direct  weighing  and 
timing  method.  Water  or  oil  was  passed  at  steady  rate  through  the  rota¬ 
meter  into  a  weighing  container • mounted  on  platform  scales.  The  time  to 
collect  a  given  weight  was  determined  by  stopwatch.  The  basic  data  are 
given  in  Table  C-l. 

For  improved  sensitivity  calibration  curves  were  drawn  in 
terms  of  mass  rate  of  flow  divided  by  rotameter  reading  versus  rotameter 
reading.  These  curves  appear  as  Figures  1-3  for,  water  and  Figures  h  -  6 


for  oil. 


ROTAMETER  CALIBRATION  DATA 


111 


Q 

C\l 

CA 

LA 

LA 

O 

p- 

ca 

CM 

-St 

O 

H 

CA 

1 

H 

H 

i — 1 

• 

O 

• 

H 

• 

O 

• 

H 

• 

CA 

NO 

o- 

CM 

H 

CM 

O 

CA 

A- 

CM 

•LA 

H 

On 

rH 

rH 

1 

rH 

0 

CA 

• 

• 

• 

• 

• 

O 

rH 

O 

1 — 1 

CA 

-S* 

CM 

O 

A- 

rH 

O 

CM 

O 

CA 

O 

CO 

rH 

O 

1 

rH 

On 

On 

• 

• 

• 

« 

• 

O 

rH 

O 

O 

CM 

•LA 

la 

O 

NO 

O 

O 

rH 

la 

On 

3 

-P 

O- 

H 

ON 

| 

O 

CO 

LA 

« 

• 

• 

• 

• 

O 

O 

O 

O 

CM 

CO 

LA 

&0Q 

AMD 

O 

O 

CM 

CA 

ON 

O 

CO 

r— 

1 

8 

On 

NO 

H 

4 

• 

• 

• 

• 

• 

• 

O 

O 

O 

0 

CM 

O 

CtJ 

© 

0) 

CO 

P3 

-P 

CM 

\ 

LA 

O 

On 

On 

CO 

©  O 

C- 

la 

NO 

A- 

LA 

H  LA 

O 

NO 

1 

O 

la 

CO 

H 

ctf 

• 

• 

• 

• 

« 

O 

O 

0 

O 

0 

rH 

CO 

TJ 

•d 

d 

O 

O 

P 

© 

s 

CA 

LA 

LA 

a* 

-P  o 

CM 

la 

NO 

NO 

•H 

Q)  -d 

0 

LA 

O 

-P 

-P 

M 

g 

• 

• 

• 

• 

• 

• 

c6 

0 

O 

CA 

O 

O 

H 

<H 

43 

O 

q 

Pi 

rH 

CO 

§ 

O 

3 

NO 

On 

CM 

S 

3 

H 

H 

CA 

0 

CA 

NO 

O 

CA 

H 

ft 

« 

• 

• 

• 

• 

• 

0 

O 

CM 

O 

O 

H 

O 

H 

-p 

CM 

1 

5 

O 

O- 

CM 

8 

a~ 

CM 

_P 

CA 

CA 

CA 

<4 

cfl 

CM 

O 

CM 

a- 

O 

CM 

A- 

© 

Ph 

• 

« 

• 

• 

• 

• 

rH 

O 

O 

rH 

O 

O 

O 

,0 

10 

ctf 

CO 

ft 3 

CA 

NO 

JEJ 

CA 

ON 

LA 

CM 

A— 

O 

d 

O 

H 

CA 

CA 

rH 

H 

-P 

•H 

HI 

O 

rH 

On 

O 

H 

CA 

• 

• 

• 

• 

• 

• 

d 

O 

O 

O 

O 

O 

O 

© 

> 

ft 

O 

•d 

LA 

H 

CA 

LA 

A- 

rH 

© 

P  ft 

O- 

A- 

A- 

C — 

c — 

A- 

d 

o*  £ 

1 

1 

1 

! 

'  1 

1 

cti 

dl 

O 

0 

O 

O 

LA 

On 

0- 

a- 

A- 

O- 

O- 

NO 

rH 

*H 

O 

•d 

u 

d 

d 

O 

1 

© 

■P 

$ 

© 

-P 

*-, 

W 

d 

$ 

$ 

d 

d 

d 

| 

• 

CM 

CM 

rH 

H 

CA 

CA 

Cti 

U  O 

C— 

LA 

A- 

LA 

A- 

LA 

O 

0)  JSi 

LA 

O 

LA 

O 

LA 

Q 

;H 

+3 

©  1 — 1 

a 

a 

CM 

rH 

_P 

rH 

a 

-St 

rH 

<H 

•H 

E  ctf 

,Q 

rO 

rQ 

£> 

& 

43 

O 

ro  *H 

LA 

CO 

rH 

CO 

LA 

CO 

© 

-P  d 

O  © 

8 

0 
c — 

O 

A- 

O 

A- 

O 

NO 

O 

A— 

% 

Pi  co 

LA 

LA 

LA 

LA 

LA 

LA 

* 

±12 


100 

90 

80 

70 

k 

60 

50 

40 

30 

20 

10 


Hj 


r 

PIG.  C. 

-2  Calibrati 
Water  Rot 
#  5708bl4 

on  of 
ameter 

052 

c 

c 

C 

> 

- 

\  p 

4  \  r 

) 

- . - Vx 

/ 

_ 

1.28  1.30  132  1.34  1.36 


w  MASS  RATE  OF  FLOW,  lb.  per  sec. 

R  ’  ROTAMETER  SETTING 


114 


FLOW  RATE ,  W ,  cu.  ft.  per  sec. 


100 

90 

80 

70 

60 

50 

40 

30 

20 

10 


115 


1.30  134  1.38 


PIG.  C-4 


Calibra 
Oil  Roi 
#  5708t 


.tion  of 
ameter 


14051 


n 


142  1.46 


Xl0+3  MASS  RATE  OF  FLOW,  lb.  per  sec. 
’  ROTAMETER  SETTING 


1.13  1.15  117  119  120 


W  |0+2  MASS  RATE  OF  FLOW,  lb.  per  sec. 
R  ’  ROTAMETER  SETTING 


117 


100 


80 


CC 


CD 


LiJ 

if) 


70 


60 


E5  50 

h- 

LU 

< 


o 

CC 


40 


30 


20 


10 


-  £ 

^ _ 

— cr 

a 

/ 

o 

PIG. 

C-6  Calibi 
Oil  Rc 
#  570? 

•at ion  of 
tame ter 
Ibl4053 

. 

V 

c 

> 

\ 

/ 

\ 

V. 

) 

p —  -0— 

i 

i 

3.50 


3.60 


370  380  3.90 

Wyi(Y2  MASS  RATE  OF  FLOW,  lb.  per  sec. 
R 


ROTAMETER  SETTING 


118 


APHSKDIX  D 


Working  Plots. 


119 


APPENDIX  D 

Appendix  D  contains  reproductions  of  the  original  "working  plots 

used  in  the  correlation  of  the  pressure  drop. 

A  typical  relationship  between  f^  and  Ry,  at  constant  D  and  Vr-.,  is 

shown  in  Figure  D-l.  This  relationship  is  basic  and  was  obtained  at 

constant  superficial  reference  phase  velocity,  Vp,  for  each  second  phase 

density,  p  .  Data  for  the  construction  of  these  relationships  is 
s 

presented  in  Appendix  B  for  second  phase  densities  of  0.092,  O.362, 

0.552  and  48.57  ibs.  per  cu.  ft«  In  addition  a  density  of  0.184  lbs. 

per  cu.  ft.  was  used  from  Short’s  (l6)  data. 

3.5 

Plots  of  fl  versus  DV-d  at  constant  D  V  ,  p  and  V  ^  are  shown  in 
in  LX  s  s  s 

Figures  D-2  and  D-6.  .  These  plots  are  constructed  from  the  basic  fp  versus 

relationships  (Figure  D-l)  in  the  following  manner.  A  value  of  the 

parameter  D''^  V  is  chosen  and  V\  is  calculated  at  the  diameter  on  the 
s  s 

fp-R  plot.  The  value  of  =  Vs  is  calculated  from  the  same  plot  since 

% 

Vp  is  constant.  Finally  the  value  of  fp  corresponding  to  Ry  is  plotted 

against  the  product  DVp.  This  procedure  is  repeated  for  all  values  of 

3.5 

VR  and  ps  at  even  increments  of  the  parameter  D"  Vg. 

2.85 

Figure  D-7  is  a  plot  of  the  product  f' (DV~)  versus  (V  ~2&)  +  a 

JK  K  S  Ja 

3.5 

constant  at  constant  Vg.  The  power  2.85  is  determined  by  considering 

•the  best  fit  slope  through  the  points  of  Figures  D-2  to  D-6.  The 
difference  in  kinematic  viscosity  (t^s  - 7/0 )  is  introduced  to  satisfy  the 
limiting  condition  for  single  phase  flow. 

In  Figure  D-8  f _L (DVp Y~ (1/  is  plotted  against  DJ*^V  . 

The  power  0.33  for  the  viscosity  is  determined  by  considering  the  range 
were  the  oil  and  air  data  overlap.  Using  the  points  in  this  region  as  a 
guide,  a  series  of  lines  sloping  at  0.33  are  constructed  through  the 
points  for  the  parameters. 
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